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Abstract We investigate the mass-to-light ratios of stellar populations as predicted by 

stellar population synthesis codes and compare those to dynamical/gravitational 
measurements. In Bell & de Jong (2001) we showed that population synthesis 
models predict a tight relation between the color and mass-to-light ratio of a stel- 
lar population. The normalization of this relation depends critically on the shape 
of the stellar IMF at the low-mass end. These faint stars contribute significantly 
to the mass, but insignificantly to the luminosity and color of a stellar system. In 
Bell & de Jong (2001) we used rotation curves to normalize the relation, but ro- 
tation curves provide only an upper limit to the stellar masses in a system. Here 
we compare stellar and dynamical masses for a range of stellar systems in order 
to constrain the mass normalization of stellar population models. We find that 
the normalization of Bell & de Jong (2001) should be lowered by about 0.05-0. 1 
dex in M/L. This is consistent with a Kroupa (2001), Chabrier (2003) or a Ken- 
nicutt (1983) IMF, but does not leave much room for other unseen components. 
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1. Introduction 

In Bell & de Jong (2001) we showed that stellar population models predict a 
strong correlation between an optical color of a stellar population and its mass- 
to-light (M/L) ratio (see also e.g., Bell et al. 2003; Portinari et al. 2004). We 
showed that the slope of this relation is rather insensitive to the exact details 
of the star formation history and chemical enrichment of the stellar popula- 
tion (except for recent star bursts) and to dust reddening, owing to the well- 



known age/metallicity/dust degeneracy. Furthermore, the color-M/L slope is 
also rather insensitive to the IMF used. Yet, the normalization of the color- 
M/L relation is highly IMF dependent, shifting up and down depending on 
how many stars are present at the low-mass end of the stellar IMF (these stars 
contribute significantly to the mass of a population, but insignificantly to its 
luminosity and color). 

In Bell & de Jong (2001) we used maximum disk rotation curves to constrain 
the normalization of the color-M/L relation. The predicted stellar population 
masses derived from the color-M/L relation should never over-predict the ob- 
served dynamical masses derived from rotation curves. However, while rotat- 
ing gas in a disk galaxy is a very simple dynamical system and hence a clean 
constraint, rotation curves of disk galaxies have the disadvantage that they only 
provide an upper limit to M/L ratios once we accept that dark matter may be 
present in disk galaxies. There is no guarantee that there is no unseen mat- 
ter contributing to the dynamical mass within the radius where the maximum 
disk is constrained, be it baryonic (e.g., cold molecular gas) or non-baryonic. 
Hence, rotation curves only provide an upper limit to the normalization of the 
color-M/L relation. Here we compare dynamical masses and masses predicted 
by stellar population modeling of a variety of stellar systems in order to con- 
strain the normalization of the color-M/L relation' . 

2. Comparing dynamical and population M/L estimates 

In order to compare dynamical and stellar population masses we have to 
make a number of assumptions: 

■ The IMFs of the stellar populations in the different objects are the same, 
notwithstanding the large range in object scale sizes and masses in- 
volved. 

■ The stellar population models used are accurate in a relative sense (not 
necessary in absolute calibration). 

■ The stellar systems in question have not selectively lost (or accreted) 
stars in a particular mass range. 

■ Where necessary we use the HST Key Project distance scale. 

We will now go through a number of dynamical/gravitational versus stellar 
population mass comparisons, and express the range of allowed population 
M/L ratios in terms of the IMF normalization used in Bell & de Jong (2001), 

' In principle, any of these systems could have a dark component co-spatial with the stellar light (in some 
cases this is rather unlikely), and hence all comparisons are strictly speaking upper limits to the relation. 
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i.e. a Salpeter x=1.35 IMF between 0.1 and 125 Mq reduced in mass by a 
factor 0.72. 

Globular Clusters At first sight globular clusters seem ideal targets to com- 
pare dynamical and stellar population masses: their stellar populations, dust 
corrections and dynamics are simple, and they are unlikely to contain large 
amounts of dark matter near their centers. However, mass segregation has re- 
sulted in centers of clusters being dominated by massive stars, the outer parts 
by lower mass stars. The outer stars are subsequently more likely to be stripped 
by interaction with the galaxy, making it even harder to get a good sampling of 
the original full IMF. Detailed dynamical modeling of Galactic globular clus- 
ters shows clear evidence of these effects, with M/L changing with radius (e.g., 
Gebhardt & Fischer 1995). 

When we compute the dynamical core M/L ratios of Galactic globular clus- 
ters following McLaughlin (2000) and compare those to single burst PEGASE 
(Fioc & Rocca-Volmerange 1997) models using the colors and metallicities 
of Harris (1996), we find that the dynamical M/L values are much lower (by 
about 0.23 dex) than predicted by the single burst models of a 10-12 Gyr old 
population. Alternatively, we can follow a more simplified approach by using 
virial masses (Pryor & Meylan 1993), which are more representative of the to- 
tal globular cluster. We find that the 12 Gyr stellar population masses are lower 
by 0.10 dex than the virial masses when using our Bell & de Jong (2001) IMF 
normalization, albeit with a large scatter of 0.20 dex rms (comparable to the 
uncertainties in the dynamical M/L values). 

In recent years it has also become possible to measure virial masses of glob- 
ular clusters in other nearby galaxies. This has the advantage that it is easier 
to get integrated properties the globular clusters and many objects at the same 
distance, but as disadvantage the limited accuracy that can be reached, even 
with 8 m class telescopes. The results of extra-galactic globular clusters are 
still inconclusive, with dynamical masses of Cen A as measured by Martini & 
Ho (2004) being more massive than our stellar population model predictions 
by 0.08 dex, but the dynamical masses of M33 (Larsen et al. 2002) being 0.27 
dex lower than predicted. 

Elliptical galaxies Recently, Cappellari et al. (2006) have performed a de- 
tailed analysis of dynamical and stellar populations masses of a sample of 
early-type galaxies. Their integral field spectrograph SAURON data allows 
them to derive accurate dynamical masses using Schwarzschild modeling and 
stellar population masses using line-strength indices modeling. Using a Kroupa 
(2001) IMF and Vazdekis et al. (1999) stellar population synthesis models, they 



^ We do not explicitly include a mass contribution for objects with masses less than 0. 1 Mq ; as argued later, 
the contribution from brown dwarf or planetary regime objects to the stellar M/L is expected to be 0.04 dex 
or less. 



find that old, fast rotating elliptical galaxies have dynamical and stellar popula- 
tion masses that are very similar. Younger, fast rotating elliptical galaxies have 
smaller stellar population masses than dynamical masses, but they can be made 
to agree by assuming that the young ages are the result of a superposition of 
a dominant, old massive population and a small, young population. However, 
slowly rotating, old massive elliptical galaxies seem to have higher dynamical 
than stellar population M/L ratios, a discrepancy that cannot be solved by a 
super-position of young and old populations, because the population akeady is 
old according to the line indices. 

Cappellari et al. (2006) argue that under the assumption that the IMF is 
the same for all galaxies this must mean that these massive, slowly rotating 
galaxies have a significant dark matter within their effective radius where the 
dynamical measure was made. However, once we accept that some elliptical 
galaxies must have a dynamically significant amount of dark matter in their 
central region, we cannot exclude that all elliptical galaxies have dark matter 
contributing to their central dynamics. Therefore, the comparison of stellar 
population and dynamical masses in elliptical galaxies becomes an upper limit 
to the normalization of the "IMF mass", identical to the maximum disk rota- 
tion curve constraint. In terms of this mass normalization, this argues for a 
~ 0.05 dex lower normalization that used by Bell & de Jong (2001), given that 
the Vazdekis models include masses down to 0.01 Mq. 

Maximum disk rotation curves: As described above, we used maximum 
disk rotation curve limits to normalize the color-M/L relation in Bell & de 
Jong (2001). In Kassin, de Jong & Weiner (2006) we have repeated this anal- 
ysis, but we expanded the Verheijen (1997) Ursa Major cluster sample with 
34 luminous galaxies and improved the treatment of shifting the mass mod- 
els to another distance. We compared the maximum disk values to the updated 
color-M/L relations of Bell et al. (2003) and find that the Bell & de Jong (2001) 
normalization is fully consistent with this expanded data set. The normaliza- 
tion may at best be 0.05 dex higher to account for the scatter in the color-M/L 
relation. 

Minimum disk rotation curves: While most galaxy rotation curves are fairly 
smooth, some show enough structure to allow determination of a lower limit 
to a stellar M/L to explain these structures under the assumption that the dark 
matter component is smooth (e.g., Noordermeer et al. 2004). Such analysis is 
complicated by the unknown intrinsic distribution of dark matter, the effect of 
adiabatic contraction, and rotation curve uncertainties (including non-circular 
motions). Using NGC 157 (Kassin et al. 2006) we find a lower limit of -0.3 dex 
with respect to the Bell & de Jong (2001) normalization to explain the strongly 
declining rotation curve of this galaxy. However, the large asymmetries and 
hence large errorbars on the rotation curve of this galaxy limits the usefulness 
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of this galaxy. More suitable systems (mainly early-type spiral galaxies with 
falling rotation curves) are studied by Noordermeer (2006). 

Disk velocity dispersions: The mass and scale height distribution determine 
the vertical velocity dispersion of a self-gravitating disk. Thus, to determine a 
disk mass from a measured vertical velocity dispersion in a face-on system we 
have to make assumptions about the (unobservable) vertical stellar distribution, 
while for edge-on systems, where we can measure the vertical stellar distribu- 
tion, we have to relate the observed radial and tangential velocity dispersion 
to the (unobservable) vertical velocity dispersion (Bottema 1997). Recently, 
Kregel et al. (2005) used velocity dispersions of a sample of 15 edge-on galax- 
ies to determine a dynamical mass TuUy-Fisher relation and compared it to the 
stellar population mass TuUy-Fisher relation of Bell & de Jong (2001). They 
find an offset of about -0.24 dex, assuming a vertical-to-radial velocity disper- 
sion ratio ipzloR) of 0.6. However, the determined M/L ratio scales quadrati- 
cally with the poorly known (Jz/(Tr ratio. To the best of our knowledge, only 
3 measurements of OzJor have been made to date, ranging between 0.5 and 
0.9 (Gerssen, Kuijken, & Merrifield 2000). Without substantially better un- 
derstanding of the behavior of ct^/o-r as a function of galaxy properties, it is 
unclear that one can place competitive constraints on stellar M/L ratios using 
this method. 

Bar streaming motions: A galactic bar moving through the interstellar medium 
creates streaming motions and often a shock, the size of which depends some- 
what on the pattern speed of the bar, but mostly on the mass of the bar. Weiner 
et al. (2001; 2004) obtained Ha velocity fields of NGC4123 and NGC3095 
and modeled these with fluid-dynamical models. Their models only permit a 
limited range in stellar M/L, such that the galaxies are close to maximum disk. 
In Fig. 1 we compare the local bar colors and the derived M/L values of these 
two galaxies to the stellar population models of a range in metallicity and with 
exponentially decaying star formation rates. We show the models normalized 
at the Bell & de Jong (2001) value on the left, reduced by 0.1 dex in M/L 
on the right. The models cover a limited area in these diagrams, showing the 
age-metallicity degeneracy that makes the color-M/L relation work in the first 
place. 

We expect the central region to suffer from extinction, and we have plot- 
ted indicative dereddening vectors on the measured data points. These vectors 
were derived from Tully et al. (1998) global galaxy reddening values, and the 
extinction in the central region may be somewhat higher. The left panel, where 
the Bell & de Jong normalization is used, shows that the raw and reddening- 
corrected stellar M/L ratios are consistent with the model normalization. The 
right-hand panel, with the model stellar M/L values decreased by 0. 1 dex com- 
pared to Bell & de Jong (2001), is just consistent with the reddening-corrected 
data. The bar streaming motions modeling therefore provides some of the 
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Figure 1. Comparing the Weiner et al. (2001; 2004) bar dynamics M/L constraints to Pegase 
stellar population models (Fioc & Rocca-Volmerange 1997). We plot 12 Gyr old exponential 
decaying star formation rate models with different decay rates r connected by solid lines, differ- 
ent metallicities by dotted lines, as indicated. The two data points for NGC 4123 and NGC 3095 
are indicated with formal errorbars. The arrows on the data points indicate the reddening for 
these galaxies according to Tully et al. (1998). On the left we show the Bell & de Jong (2001) 
normalization, at the right a 0.1 dex lower normalization. 



Strongest constraints on our color-M/L normalization, allowing only a range 
of ~0.2 dex. 

Perez et al. (2004) confirm the analysis of Weiner et al. (2001) to the extent 
that, for the two out of their sample of five galaxies for which they could derive 
M/L constraints, the barred galaxies had to be close to maximum disk (at least 
80% stellar mass contribution in bar region). 

Spiral arm streaming motions: In a similar fashion, streaming motions can 
be used to estimate the mass in a spiral density wave. Clearly this is a more 
challenging exercise, as arm-induced shocks and streaming motions are much 
weaker than those induced by bars. Kranz et al. (2003) studied five high surface 
brightness galaxies with long-slit spectra and optical/near-IR surface photom- 
etry. They can only weakly constrain stellar M/L (their Table 3), and find /) 
most of their sample have maximum disk M/L values consistent with the Bell 
& de Jong calibration, and ii) most high Urot disks are consistent with close 
to maximum disk, whereas the lower rotation velocity disks could be substan- 
tially sub-maximal (with less than ~60% of the disk mass coming from stars 
within 2.2 disk scalelengths). 

Strong galaxy lensing: Strong gravitational lensing provides an estimate of 
the gravitational mass within the lens area with rather straightforward model- 
ing. To compare the gravitational mass of the lens with its stellar population 
mass we have to correct the observed colors using k-corrections or better yet, 
redshift the model spectra to the lens redshift and calculate a new color-M/L 
grid in the observed bands. Furthermore, we have to realize that galaxies are 
younger at higher redshift and the color-M/L relation will shift. For simplic- 
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ity we can model this with exponentially declining star formation rate models 
(which becomes an increasingly poorer approximation at higher redshift, be- 
cause starbursts will become relatively more important). Such exponentially 
declining models — started 12Gyr ago — have color-M/L relations in rest- 
frame colors that are decreased by 0.15 dex at redshift 1 compared to the z=0 
relation. 

To minimize these corrections a lensing galaxy at low redshift should be 
used in the gravitational to stellar mass comparison. Furthermore, the lensed 
images should be of small angular separation, surrounding only the central 
region of the lens galaxy, which is most likely to be dominated by stellar mass. 
We should keep in mind that any masses derived from lensing, like maximum 
rotation curves, only provides an upper limit to the stellar population mass 
estimates, as significant dark matter may be present in the centers of some 
galaxies as also suggested by the Cappellari et al. (2006) elliptical galaxies 
result. 

Indications from the first studies to satisfy these criteria are encouraging. 
Smith et al. (2005) present an example of a low redshift (^=0.0345), tight lens, 
finding a M/L/~1.8, M/Lb~4.7. Using the reported F475W and F814W 
observed magnitudes, galactic foreground corrections, and k-corrections as- 
suming a non-evolving ancient galaxy template (solar metallicity and ~12 Gyr 
old), we find excellent agreement with the predicted stellar M/L values, as usual 
accounting for gas recycling from ageing stellar populations (M/L/ pi.cd~l-8» 
M/Lg pi.ed'^4.5). Koopmans et al. (2006) analyze an extensive sample of 15 
lenses with z=0.06-0.33. In this case, the dynamically-derived total M/L scale 
is compared with the lensing results, finding consistency (i.e., they have roughly 
cross-checked the dynamical mass scale — e.g., Cappellari et al.'s scale — 
with the lensing scale). A more careful, explicit test of the color-derived stellar 
mass scale with the lensing mass scale is clearly warranted. 

3. Conclusions 

In Fig. 2 we give an overview of all constraints derived in the previous sec- 
tion on stellar population M/L values relative to an IMF normalization of Bell 
& de Jong (2001). The strongest constraints are currently provided by the 
Weiner et al. (2001, 2004) constraints from bar streaming motions. However, 
this constraint is derived from only two galaxies, and is therefore very suscep- 
tible to for instance errors in the distances to the galaxies. We have indicated 
in Fig. 2 the effect of 15% distance errors, which accounts for 0.075 dex off- 
set in the M/L IMF normalization. Another source of uncertainty lies in the 
stellar population modeling. The Padova isochrone tracks (Girardi et al. 2002) 
are used in many of the popular models (e.g.. Chariot & Bruzual 2003; Fioc 
& Rocca-Volmerange 1997; Vazdekis 1999) and these models all give very 



1 1 1 1 1 


1 1 1 1 > 1 1 


error 
Lensing 
Spiral arms 
Bar stream motions 
Disk vel. dispersion 
Min disk Rot Curve 
Max disk Rot Curve 
Eliptical Gals 
Extra-Gal GC 
Galactic GC 

1 1 1 1 1 1 1 






































1 1 1 1 1 


1 1 1 1 1 1 1 



0.20 0.00 -0.20 -0.40 -0.60 

Alog(M/L) (Bell 8c de Jong) 



-0.80 



Figure 2. Overview of offsets between dynamical and stellar population mass estimates rela- 
tive to the normalization of Bell & de Jong. Ranges are indicated by horizontal errorbars, upper 
and lower limits are indicated with horizontal arrows. The vertical dotted line indicates the 
our best estimate for the offset given all constraints. The uncertainty in this offset caused by 
15% distance errors are indicated at the top. The vertical arrows indicate offsets expected for a 
Kroupa (2001) IMF (left, green) and a Kennicutt (1983) IMF (right, cyan). 



similar results in terms of the color-M/L relation. However, a class of mod- 
els based on the fuel consumption theorem that treats the late stages of stellar 
evolution differently are giving somewhat different results, especially in the 
infrared where AGB and RGB stars dominate (e.g., Maraston 2005). These 
models have indeed a slightly steeper slope in the color-M/L relation, espe- 
cially in the near-infraied, but surprisingly the normalizations of the different 
sets of models are very consistent in the intermediate color range where most 
of the normalization constraints are derived. Combining all constraints, we 
argue that the most likely normalization is about 0.05-0.1 dex lower than the 
Bell & de Jong (2001) normalization, i.e. a Salpeter x=1.35 IMF between 0.1 
and 125 Mq reduced in mass by a factor 0.6 (-0.22 dex). 

This normalization is consistent with for instance the Kroupa (2001) or 
Chabrier (2003) IMFs (offset by about -0.25 dex with respect to Salpeter IMF) 
and a Kennicutt (1983) IMF (-0.3 dex). It leaves however not much margin 
for other unseen but known mass components in the dynamical comparisons. 
Stellar and substellar objects below our 0.1 Mq mass limit could contribute up 
to 10% (0.04 dex) of the IMF mass. A significant molecular gas component in 
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the inner part of spiral galaxies could push down the rotation curve constraint. 
However, it is satisfying to see that the current Galactic IMF estimates give 
stellar mass estimates in a wide range of objects that are fully consistent with 
their dynamical masses. 

References 

Bell, E. R, & de Jong, R. S. 2001, ApJ, 550, 212 

Bell, E. P., Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289 

Bottema, R. 1997, A&A, 328, 517 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 

Cappellari, M., et al. 2006, MNRAS, 366, 1 126 

Chabrier, G. 2003, ApJ, 586, L133 

Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 

Gebhardt, K., & Fischer, R 1995, AJ, 109, 209 

Gerssen, J., Kuijken, K., & Merrifield, M. R. 2000, MNRAS, 317, 545 

Girardi, L., Bertelli, G., Bressan, A., Chiosi, C., Groenewegen, M. A. T, Marigo, R, Salasnich, 

B., & Weiss, A. 2002, A&A, 391, 195 
Harris, W.E. 1996, AJ, 112 

Kassin, S. A., de Jong, R. S., & Weiner, B. J. 2006, ApJ in press (astro-ph/0602027) 
Kennicutt, R. C, Jr. 1983, ApJ, 272, 54 
Koopmans, L.V.E., Treu, T., Bolton, A.S., Buries, S., Moustakas, L.A., 2006, submitted to ApJ 

(astro-ph/0601628) 
Kranz, T., Slyz, A., & Rix, H.-W. 2003, ApJ, 586, 143 
Kregel, M., van der Kruit, R C, & Freeman, K. C. 2005, MNRAS, 358, 503 
Kroupa, R 2001, MNRAS, 322, 231 

Larsen, S. S., Brodie, J. R, Sarajedini, A., & Huchra, J. R 2002, AJ, 124, 2615 
Maraston, C. 2005, MNRAS, 362, 799 
Martini, R, & Ho, L. C. 2004, ApJ 610, 233 
McLaughlin 2000, ApJ, 539, 618 

Noordermeer, E. 2006, Ph.D. Thesis, University of Groningen 

Noordermeer, E., van der Hulst, T, Sancisi, R., & Swaters, R. 2004, lAU Symposium, 220, 287 
Perez, I., Fux, R., & Freeman, K. 2004, A&A, 424, 799 
Portinari, L., Sommer-Larsen, J., & Tantalo, R. 2004, MNRAS, 347, 691 
Pryor, C., & Meylan, G. 1993, in ASP Conf. Ser. 50, Structure and Dynamics of Globular 

Clusters, ed. S. Djorgovski & G. Meylan (San Francisco: ASP), 357 
Salpeter, E. E. 1955, ApJ, 121, 161 

Smith, R. J., Blakeslee, J. R, Lucey, J. R., & Tonry, J. 2005, ApJL, 625, L103 
Tully, R. B., Pierce, M. J., Huang, J.-S., Saunders, W., Verheijen, M. A. W., & Witchalls, R L. 

1998, AJ, 115,2264 
Vazdekis A., 1999, ApJ, 513, 224 

Verheijen, M. A. W. 1997, Ph.D. Thesis, University of groningen 
Weiner, B. J. 2004, lAU Symposium, 220, 265 
Weiner, B. J., Sellwood, J. A., & Williams, T. B. 2001, ApJ, 546, 931 



